LBNL-56660 



Multiple Parton Scattering in Nuclei: Heavy Quark Energy Loss and 
Modified Fragmentation Functions 

Ben- Wei Zhang^, Enke Wang^^^ and Xin-Nian Wang^ 
^Institute of Particle Physics, Huazhong Normal University, Wuhan 430079, China 
^Nuclear Science Division, MS 70R0319, Lawrence Berkeley National Laboratory, Berkeley, CA §4720 USA 

Multiple scattering, induced radiative energy loss and modified fragmentation 
' functions of a heavy quark in nuclear matter are studied within the framework of 

generalized factorization in perturbative QCD. Modified heavy quark fragmentation 
functions and energy loss are derived in detail with illustration of the mass dependen- 
cies of the Landau-Pomeranchuk-Migdal interference effects and heavy quark energy 
loss. Due to the quark mass dependence of the gluon formation time, the nuclear 
size dependencies of nuclear modification of the heavy quark fragmentation function 
and heavy quark energy loss are found to change from a linear to a quadratic form 
when the initial energy and momentum scale are increased relative to the quark mass. 
The radiative energy loss of the heavy quark is also significantly suppressed due to 
^ . limited cone of gluon radiation imposed by the mass. Medium modification of the 

I heavy quark fragmentation functions is found to be limited to the large z region due 

. to the form of heavy quark fragmentation functions in vacuum. 

o 



O 

o 



I. INTRODUCTION 



o 

Ph In ultra-relativistic heavy-ion collision, energetic partons from initial hard processes have to prop- 

CD . 

(-H ^ agate through the produced dense medium and therefore suffer multiple scattering and lose a sig- 

nificant amount of energy. Such energy loss or jet quenching has been proposed as a good probe of 
I the hot and dense medium formed in high-energy nuclear collisions [1,2]. It in effect suppresses the 

^ I final leading hadron distribution from the propagating parton giving rise to modified fragmentation 

functions and the final hadron spectra [3,4]. Recent theoretical studies [5-9] all show that the effec- 
tive parton energy loss is proportional to the gluon density of the medium. Therefore measurements 
of the parton energy loss will enable one to extract the initial gluon density of the produced hot 
medium. Strong suppression of high transverse momentum hadron spectra is indeed observed by 
experiments [10-13] at the Relativistic Heavy-Ion Collider (RHIC) at the Brookhaven National Lab- 
oratory (BNL). The suppression pattern agrees very well with the jet quenching mechanism [14,15], 
indicating large parton energy loss in a medium with large initial gluon density. Comparing to jet 
quenching as measured in deeply inelastic scattering off nuclei, the initial gluon density in central 
Au + Au collisions at ^/s = 200 GeV is about 30 times higher than that in a cold Au nuclei [16]. 
Such a high initial gluon density is unprecedented and is a strong indication of the formation of 
quark gluon plasma. 

The extraction of the initial gluon density from jet quenching pattern as measured at RHIC relies 
on the assumption that it is caused by parton multiple scattering and induced radiation. While such 
an assumption is based on a solid physical picture and is supported by a multitude of experimental 
data [17], it is still important to have additional and independent study of the consequences of parton 
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energy loss. Quenching of heavy quark spectra has been proposed as a special probe because of the 
unique mass dependence of the energy loss and medium modification of the fragmentation functions 
[18-21]. In this paper, we will apply the framework of twist expansion that was developed for the 
study of medium modification of parton fragmentation functions to heavy quarks. 

Multiple parton scattering inside nuclei in generally is a higher twist process that involves multiple 
parton correlations. By dimensional counting, such higher twist processes are power-suppressed in 
terms of the momentum scale Q"^ involved. Therefore, one can have a systematic expansion of the 
cross section in l/Q^. Since the probability of multiple scattering increases with the nuclear size, 
the leading higher twist contribution should be enhanced by A}-!'^ . Furthermore, multiple parton 
correlations involve only the intrinsic properties of the nuclei. They should be independent of the 
hard processes involved. This twist expansion is referred to as generalized factorization [22]. Such 
framework has been applied to semi-inclusive DIS on nuclear targets to study nuclear modification 
of light (massless) quark fragmentation functions and effective parton energy loss [23]. Because 
of the Landau-Pomcranchuk-Migdal interference [24], the phase space available for the induced 
gluon radiation is limited that is also proportional to The final twist-four contribution to the 

modified fragmentation functions is then proportional to A^/^/Q^ [23]. Such a quadratic dependence 
on nuclear size is indeed observed in semi-inclusive deep inelastic lepton-nucleus experiments [16,25]. 

In this paper, we will extend the study of modified fragmentation functions in nuclei [23,26] to 
heavy quarks when they propagate through nuclear matter. We will derive the modified heavy quark 
fragmentation functions and the effective energy loss. To demonstrate the effect of quark mass, we 
will compare the results with the ones for light quarks. One of the most important effect is the 
reduction of gluon formation time when it is radiated from a slow heavy quark whose virtuality is 
not much larger than its mass. Such a reduction will effectively eliminate LPM effect and the nuclear 
size dependence of the modification will become linear in contrast to the case of a light quark. The 
second mass effect is the significant reduction of induced quark energy loss due to limited gluon 
radiation angle imposed by the mass. With detailed data analysis of experimental data both in 
the single electron channel [27] and direct measurement of heavy mesons [28], one could learn more 
about the parton energy loss mechanism in dense matter. 

The results of heavy quark energy loss in the present study were already reported in Ref. [20] . In 
this work, we will elaborate the detailed derivation and focus on numerical calculations and discus- 
sions about modified heavy quark fragmentation functions. The paper is organized as follows. In the 
next section we will present the theoretical formalism of our calculation including the generalized 
factorization of twist-4 processes. In Section III we will show in detail the calculation of different 
contributions to the modified heavy quark fragmentation function and energy loss due to multiple 
scattering. In Section IV we will numerically evaluate and discuss the modified fragmentation func- 
tions of a heavy quark propagating in nuclei. In Section V, discussion and numerical calculation of 
heavy quark energy loss will be given. We will demonstrate the mass effects by discussing how the 
dependence on medium size changes from a linear to a quadratic dependence when the energy of 
the heavy quark and the momentum scale is increased, and the suppression of the energy loss for 
the heavy quark relative to a light quark due to "dead-cone" effect [18]. Section VI will summarize 
our work. 
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II. GENERALIZED FACTORIZATION 



In order to separate the complication of heavy quark production and propagation, we consider 
a simple process of charm quark production via the charge-current interaction in DIS off a large 
nucleus. The results can be easily extended to heavy quark propagation in a hot medium. The 
differential cross section for the semi-inclusive process i{Li) + A{p) — > ve.{L2) + H{£h) + X can be 
expressed as 

Here Li and L2 are the four momenta of the incoming lepton and the outgoing neutrino, £h the 
observed heavy meson momentum, p = m^/2p+, 0_l] is the momentum per nucleon in the 
nucleus, mjv is the mass of nucleon and s = {p + ii)^. Gp is the four-fermion coupling constant, 
and q = L2 — Li = [—Q'^/2q~,q~,0±] the momentum transfer via the exchange of a vector boson 
B{q). The charge-current leptonic tensor is given by = l/2Tr(/,i7^(l — 75) ,£2(1 + 75)7i^)- We 
assume <C M^. The semi-inclusive hadronic tensor is defined as, 

^ = I H){X, HU+^W-Kd^q+p- px - Ih) (2) 



X 



where runs over all possible final states and J+ = '^'fluVipf is the hadronic charged current. 
Here we define 1^ = (1 — 75)!^^ and Vij represents the CKM flavor mixing matrix [29]. We want to 
clarify that the symbol Q in this paper stands for both the heavy quark flavor and the momentum 
scale of the exchanged vector boson. 

Wc consider the process of DIS in which collides with a light quark q with momentum k — xp 
producing a heavy quark Q with mass M and momentum £q. The heavy quark then fragments into 
a heavy quark flavored hadron. In this paper, we will take the charm quark as an example and the 
cases for other heavy quarks will be straightforward. In order to investigate the energy spectrum of 
charm fragmentation we define the Lorentz-invariants z = iQ/q~ ,zh = ^h/^c- The leading-twist 
and lowest order perturbative QCD (pQCD) calculation of heavy quark production gives 



dW^ 

dZH 



= Y,j dxf^{x)Hj°){x,p,q,M)DQ^HizH) , (3) 



where f^{x) is the quark distribution and Dq^h{zh) is the non-perturbative heavy quark frag- 
mentation function in vacuum [30-33]. The hard partonic part is 

q,M) = ^ Tr{x]h^.V{fi+x]^)Vh,) ^S{x - xb - xm) , (4) 

In the case when the momentum scale Q is much larger than the heavy quark mass, large loga- 
rithms such as log(Q^/M^) arise to all orders of the perturbative expansion, so the fix-order perturba- 
tion theory breaks down and a perturbative resummation of large quasi-collinear logs, log((3^/M^), 
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should be performed [31,32], which wiU give the corresponding QCD evolution equations for par- 
ton distribution functions and heavy quark fragmentation functions. After considering higher order 
contributions, the inclusive tensor can be written as [31,34-37] 

-J^^Y. J dxf^{x,fi')Hj^J{x,p,q,M)DQ^HizH,t^'), (6) 

where DQ^if{zff,fi^) satisfies the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [38] QCD 
evolution equations. 

In a nuclear medium, the propagating heavy quark in DIS will experience additional scattering with 
other partons from the nucleus. The rescattering may induce additional gluon radiation and cause the 
heavy quark to lose energy. Such induced gluon radiation will effectively give rise to additional terms 
in the evolution equation leading to the modification of the heavy quark fragmentation functions in 
a medium. These are higher-twist corrections since they involve higher-twist parton matrix elements 
and are power-suppressed. We will consider those contributions that involve two-parton correlations 
from two different nucleons inside the nucleus. They are proportional to the size of the nucleus [39] 
and thus are enhanced by a nuclear factor A^^^ as compared to two-parton correlations in a nucleon. 
As in previous studies [23] , we will neglect those contributions that are not enhanced by the nuclear 
medium. 




FIG. 1. A typical diagram for quark-gluon re-scattering processes with three possible cuts, central(C), 
left(L) and right(R). 



For the production of heavy quarks there are usually two kinds of mechanisms: intrinsic and 
extrinsic (via gluon fusion) heavy quark production [40] . Since we are only interested in rescattering 
and induced gluon radiation of the heavy quark after it is produced in DIS we will only consider 
a simple case of intrinsic heavy quark production via charged-current weak interaction. In this 
case, after small modification we can still employ the generalized factorization of multiple scattering 
processes [22] . In this formalism with collinear approximation, the double scattering contribution to 
radiative correction from processes like the one illustrated in Fig. 1 can be written in the following 
form. 



dW, 



D 



dzH 



E 
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dz 



Dq^h{zh/z) 



dy 
27r 



'^Vi T^d^^rH^i^iv ,yi ,y2 ,kT,p,q,M,z) 



(27r)2 



it„.ir„ 1 



X e-^'^^-y-- {A\^g{0) 7+ A+{y^,OT) A+{y^,yT) Mv'M) 



(7) 
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Here H^^[y ,y-^ ,kT,p,q,M,z) is the Fourier transform of the partonic hard part 
Hij,v{x, Xi,X2, kT,p, q, M, z) in momentum space, 

xH^^{x,xi,X2,kT,p,q,M,z), (8) 

where k = [xp'^, 0, Oj_], fci = [xip'^, 0, 0_l], ^2 = [x2P'^, 0, kr], ks = [x^p'^, 0, kr], and kr is the relative 
transverse momentum carried by the second parton in the double scattering. We assimie that kr is 
small and therefore can make collinear expansion of the hard partonic cross section with respect to 
the transverse momentum of the initial partons. The first term in the collinear expansion gives the 
eikonal contribution to the leading-twist results, making the matrix clement in the single scattering 
process gauge invariant, while the second (or linear) term vanishes for unpolarized initial and final 
states after integration over kr- The leading term in the collinear expansion that contributes to the 
double scattering process comes from the quadratic term in the collinear expansion, 

^ = E/ ^Dq^h{zh/z) J ^dy^dy^l{A\i^,{0)j+F„+{y^)F+'^{y^)%{y-)\A) 



dZH 

/ 1 A r 1 

(9) 



1 92 -y^D , _ _ _ , ,^ , 

2 i^"^'^ ,yi ,y2 ^kT,p,q,M,z) 



There are many diagrams involving double parton scattering. The hard part of the partonic 
scattering for each diagram, H^^{x, xi,X2, kT,p, q, M, z), always contains two ^-functions from the 
on-shell conditions of the two cut-propagators. These 5-functions, together with the contour inte- 
grations which contain different sets of poles in the un-cut propagators, will determine the values 
of the momentum fractions x,xi, and X2 [23]. The phase factors in H^j^{y~ ,y^ ,y2 ,kT,p,q, M, z) 
[Eq. (8)] can then be factored out, which will be combined with the partonic fields in Eq. (9) to form 
twist-four partonic matrix elements or two-parton correlations. The double scattering corrections 
in Eq. (9) can then be factorizcd into the product of fragmentation functions, twist-four partonic 
matrix elements and the partonic hard scattering cross section. 

III. DOUBLE SCATTERING AND INDUCED GLUON RADIATION OFF A HEAVY 

QUARK 

According to the generalized factorization theorem in Eq. (9), we should calculate the hard part of 
parton multiple scattering. We will assume an axial gauge n ■ A — with n — [1,0^, 0^]. The hard 
part of quark-gluon double scattering has a total of 23 cut diagrams as illustrated in Figs. 11-21. 

We take the central cut diagram in Fig. 1 as an example to show how to calculate the hard part. 
With the conventional Feynman rule in the Standard Model one can write down the hard partonic 
part of the central cut-diagram of Fig. 1, 

Hc^Ay-,yi,y2,kT,p,q,z) = J d^^^e"^f^''-+'-=^^«r+^(--— y ^ 

^x^hf^Vn'^nPH^pV^'y^ 2n6+{f) 6{1 - z - ^) . (10) 



X - — Tr 
2 



q~ 
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Cf 4 7-(g + fei) + M 7 ■ (g + fci - ^) + M , ^ 



2Af/ (g + /ci)2-M2 



{q + h- If -M'^-ic 



(11) 



{q + k-£f-M'^+ ie [q + kf 

where = = 0, p"^ = m%, £q = M^, ttin and M are the nucleon mass and the produced heavy 
quark mass, respectively. In addition, is the polarization tensor of a gluon propagator in 

the axial gauge and i, £q = q + ki + k2 — i are the 4-momenta carried by the gluon and the final 
heavy quark, respectively, z = (■Q/q~ is the fraction of longitudinal momentum (the large minus 
component) carried by the final heavy quark after gluon radiation. 

In order to simplify the calculation of the trace part and extract the leading contribution in the 
limit £t ^ and /cy ^ 0, we also apply the collinear approximation to complete the trace of the 
product of 7-matrices as 



(12) 



Therefore we have 



1„ , , 4. n 1 



Tr [xin^Vin ■ + M)V^^] — Tr j-n'^H^pU" 



(13) 



After carrying out momentum integration in x, Xi, X2 and £^ with the help of contour integration 
and ^-functions, the partonic hard part can be factorized into the product of the matrix Hj^J{k,q) 
of vector boson and quark scattering, and the quark-gluon rescattering part , 

H^,,iy^,yi,y2^kT,p,q,M,z) = J dxHfJ{k,q,M) h'^ {y~ ,kT,x,p,q,M,z) , (14) 

where HjS) {k,q, M) is defined in Eq. (4). Contributions from all the diagrams have this factorized 
from. Therefore, we need only list the rescattering part for different diagrams in the following. 
We also define the momentum fractions 



XL = 



, xn = 



k^ - 2kT ■ ix 



2p+q~z{l — z) 2p+q~z 
as in the light quark case. Then, for the central-cut diagram in Fig. 1 we have [23], 

(1 + Z2)4 + (1 - z)4M2 



(15) 



, j/2 , kx, x,p, q, M,z) = J 



z){i^ + {l-z)^M^y 



X -^Cf ^;^Ii^c{y~ ,yi ,y2 JT,kT,x,p,q,M,z) , (16) 



h,c{y- ,yr ,y2 AT,kT,x,p,q,M,z) = e^(-+--)p+''"+-°f+(''i"-f2-)^(-y2-)%- - yf) 

I _ g-i(xL + {l-z)xM/z)p+y~ _ ^-i{xL + {i-z)xM/z)p+{y~-y~) 



(17) 



The above contribution resembles the cross section of dipolc scattering and contains essentially 
four terms. The first diagonal term corresponds to the so-called hard-soft process where the gluon 
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radiation is induced by the hard scattering between the vector boson B and an initial quark with 
momentum k. The quark is knocked off-shell by the B boson and becomes on-shell again after 

radiating a gluon. Afterwards the on-shell quark (or the radiated gluon) will have a secondary 
scattering with another soft gluon from the nucleus. The second diagonal term is due to the so- 
called double hard process where the quark is on-shell after the first hard scattering with the vector 
boson. The gluon radiation is then induced by the scattering of the quark with another gluon that 
carries finite momentum fraction xl + {^ — z)xm / z+xd- The two off-diagonal terms are interferences 
between hard-soft and double hard processes. The cancellation between the two diagonal and off- 
diagonal terms essentially gives rise to the LPM interference for the induced gluon radiation. Prom 
Eq. (A2), we get the formation time for gluon radiation from a heavy quark, 

= {XL + (1 - z)xmIz)p+ ' ^^^^ 

In the limit of collinear radiation (xl ^ 0) or when the formation time of the gluon radiation, r^, 
is much larger than the nuclear size, the two processes (soft-hard and double hard) have destructive 
interference, leading to the LPM interference effect [24]. It is interesting to note that the formation 
time of gluon radiation off a heavy quark is shorter than that off a light quark [23] 

rj^^. (19) 

^ XLP+ 

This is simply because the formation time for gluon radiation is only relative to the propagation of 
the quark. The velocity of low energy heavy quarks is much smaller than that of a light quark. The 
corresponding gluon formation time is also smaller. One should then expect the LPM interference 
effect to be significantly reduced for intermediate energy heavy quarks. It can be shown that this 
phenomenon will also modify the dependence of the heavy quark energy loss on the nuclear size, 
which will be discussed in detail in Section V. 

In addition to the central-cut diagram, we also should take into account the asymmetrical cut- 
diagrams in Fig. 1, which represent interference between single and triple scattering. We note that 
the trace part is the same as in the central-cut diagram. The difference is just the phase factor. 
Thus we have from these cut diagrams, 

h,L{y~,yi,y2JT,kT,x,p,q,M,z) = -e^i-+--)p^v^+^--p^(y^-y^)e{y^ -y^)e{y- -y^) 

X (^l — Q-ii^L + (l-z)xM/z)p'^{v^ -Vi)^ ^ ^20) 

hAy~,yi,y2JT,kT,x,p,q,M,z) = -e'(-+--)p^y-+^-opHv-,-v^)e{-y-)e{y- -y-) 

(^l _ g-i{xL + {l-z)xM/z)p+y~ ^ ^ (21) 

With the same procedure we can obtain the contributions of all other cut diagrams, which are 
listed in Appendix A. 

IV. MODIFIED HEAVY QUARK FRAGMENTATION FUNCTION 

To calculate the leading twist-four contribution to the semi-inclusive cross section according to 
the generalized factorization formula in Eq. (7), one has to expand the hard partonic cross section 
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in kr- We can rearrange different contributions according to the ^-functions in ^ (the sum of 
all the contributions from central-cut, right-cut or left-cut diagrams) and define 



H 



X [HS9{-y2)e{y- - yf) - H^e{-y2)9{y2 -y^)- HEO{y- - y^)0{y^ - y^)] 
^Prom Eqs. (A1)-(A36), we can check that 



HcikT = 0) = Hj,{kT =0)=H^ {kr = Q)Cf 



(1 + z")4 + (1-z)4M^ 



(22) 



(23) 



Since one can reorganize the ^-functions as 

dy^dy^ [e{-y2)0{y2 - Vi) + 0{y- - y^)e{y^ - y^) - e{-y2)0{y~ - Vi)] 

rv~ rVi 
= I dy^ I dy2 , 
Jo Jo 



I- 



(24) 



one finds that H [kx = 0) in Eq. (7) just gives the eikonal contribution to the next-leading-order 
correction of single scattering which can be gauged away since it does not correspond to any physical 
double scattering. 

The leading contributions to the quark-gluon rescattering result from the quadratic term in the 
kr expansion of , 



^lT^C{L,R)\kT=0 = ^Ca — 



1 + z' 



HDc{L,R)+OiXB/Q'e'T) 



(25) 



z [ll + (1 - zfM^]^ 

JjD _ ci[z,i^,M'^){l — e^*(^'^ + (^~^)'^"/^)P'''^2')(l — g^-^i^L + (^-z)xM / z)p+ (y- -y~)-^ 
+ C2{z, M^) Q-i{xL + {l-z)xM/z)p+y:^ ^-^ _ ^-i(xL + (l-z)xM / z)p+ {y^ -y^)^^ 

+ cz{z,i%, M^)e~'^^^"'"^"^~^^'^'^/^^^^^^~~^i~^e~'^^^"'"^^~^^^'^/^^^^''2' (26) 

JjL _ ^^^_j,^£2 ^^2^^g-i(xi, + (l-2)xM/2)p+(!/~-J/r) _ g^-i(xL + ('i.-z)xM / z)p-^ (y- -y~)-^ 

+ Cr,{z,i%,M'^){l - e-^^^L+{^-^>M/z)p+{y--y-)-^ (27) 

JjR _ c4(2;,£|.,M^)(e~'^^^ + '^~^^^'^/^^^^''2' _ Q-i(xL + ('i^-z)xM/z)p+y~)^ 
+ C5(0,4,M2)(1 - e-'ixL+0-z)xM/z)p+y-^^^ 



where the coefficient Ci{z, M^), i = 1,2, 3, 4, 5 are polynomial functions of jv^ 



cr{z,l\,M-') 

C2{Z,4,M^) = 

cs{z,£%,M^) = 



{I- zf{z'^ -Qz + l)M'^ 2z{l-zfM'^ 



0-2 



1 + Z2 4 ' 



(1-z) 



1 + Z2 

{I- z) {2z^ - 5z + 8z - 1) 2Cf 



(1 + ^2) 



+ 



Ca 



{l-zf 



z{l- z)^{3z-l) 2CFil-zy 



(l + z2) 



Ca (l + z^) 



M4 



Cpjl - zf 
Ca 



1 - 



8z(l - zf M2 (1 - zf{z^ - 4z + 1) 



1 -h ^2 £2 



1 + Z^ 



(28) 
(29) 

(30) 
(31) 
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) = {l-z)- 



(1- 



M2 



1 + 



1 

(1 



Ca 



(1 -zf + 2z-l 



zf M2 



1 + ^2 



(32) 
(33) 



In the limit of = 0, one can recover from Eqs. (25)-(28) the results for light quark multiple 
scattering with complete calculation beyond helicity amplitude approximation [26] . 

Compared with the results of gluon radiation from the light quark multiple scattering [23,26] 
which have an overall form l/^f^, the radiative gluon spectrum from a heavy quark in Eq. (25) is 
suppressed by a factor 



/q/9 



1% + (1 - 0)2M2 



1 ^ 









(34) 



where = M/q~ and 6 = ItII~ is the angle of the radiated gluon relative to the heavy quark. 
One can see that the mass of the heavy quark provides a lower bound for the radiation angle of the 
collinear gluon which dominates the gluon spectrum in the case of radiation off a light quark. This 
effectively suppresses the gluon radiation at angle smaller than the ratio of the quark mass M to its 
energy q~ , and thus reduces radiative energy loss of a heavy quark. Such a suppression of small angle 
gluons is often referred to as the "dead cone" effect [18]. In the result of our current calculation, 
additional mass effects on the final gluon spectrum arc in the mass dependence of coefficient functions 
Ci(z, M^), i = 1, 2, 3, 4, 5 and most importantly the mass dependence of the gluon formation time 
[Eq. (18)] which will dictate the LPM interference pattern in induced gluon bremsstrahlung off 
a heavy quark. These additional mass effects will significantly influence the final radiative gluon 
spectra and modified the so-called "dead cone" effect. In fact, the coefficient functions Ci{z,£'^, M'^) 
contain terms like t. Multiplied with the factor in Eq. (34), they give finite contribution to the 
gluon spectra at It = 0. This in effect will fill up the small angle cone along the direction of the 
propagating heavy quark with soft gluons, as also pointed out in Ref. [21]. The net "dead cone" 
effect still results in the significant reduction of the heavy quark energy loss, as numerical results 
shown in Fig. 9 and Fig. 10 of Section V below. 

Substituting Eqs. (25)- (27) in Eq. (9), we have the leading higher-twist contribution to the semi- 
inclusive tensor of heavy quark fragmentation in DIS off a nucleus. 



2'Ka, 



X ——^T^''^{x,XL,M^) + {g- fragmentation) 



[4 + (1 - Z)2M2]4 ^ 

(virtual corrections) , 



where 



T^f{x,XL,M^) = I ^-dy^dy^Hg\{A\^,{Q)^+F+{y^)F+%y^)^p,{y-)\A) 
T^^'\x,XL,M^) = I ^dy^dy^HEl{A\^lj,{0)'y+F,+ {y^)F+^y^)My-)\A) 



(35) 



(36) 



(37) 
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X e'(-+--)P+^"^(2/- - yr)^(yr - y^) , (38) 

X e'^-+^-^-P*y- e{-y^)e{y^ - yr) , (39) 

are twist-four quark-gluon correction functions inside the nucleus with H§ , H£ and given in 
Eqs. (26)-(28). They are all independent with each other because they involve different functions. 
The twist-four parton matrices Tj^'^ and T^g^ are new, which purely result from the mass effect 
of the heavy quark and involve left and right cut diagrams. These new parton matrices will vanish 
when we take = [See the definitions in Eqs. (26-28)]. Furthermore, they are proportional to 
coefficients 04(2, £^,M^) and 05(2, ^|,,M^) which in turn contain an additional factor of (1 — 
In the case of soft gluon radiation z ^ 1, these matrix elements are suppressed compared to the 
matrix elements involved in the central cut diagrams. 

During the coUinear expansion, we have kept (.t finite and took the limit Ut — > 0. As a conse- 
quence, the gluon field in one of the twist-four parton matrix elements in Eqs. (37)-(39) carries zero 
momentum in the soft-hard process. As argued in Ref. [23] , this is due to the omission of higher order 
terms in the coUinear expansion. As a remedy to the problem, a subset of the higher-twist terms in 
the coUinear expansion can be resummcd to restore the phase factors such as exp(ia;Tj3^y ), where 
Xt = {k^)/2p'^q~ z is related to the intrinsic transverse momentum of the initial partons. Therefore, 
soft gluon fields in the parton matrix elements will carry a fractional momentum xt- 

Until now we have only considered quark-gluon double scattering in a nucleus. To make a complete 
calculation we should also take into account the processes of quark-quark scattering. However, it 
has been shown [23] that the contributions of quark-quark scattering are suppressed by a factor 
l/Q^ as compared with quark-gluon double scattering. In the heavy quark case, it further involves 
intrinsic heavy quark distribution inside the nucleus which should be very small as compared to light 
quark and gluon distributions. Thus, we can completely neglect the contributions of quark-quark 
scattering for heavy quark propagation. 

The virtual corrections in Eq. (35) can be obtained via unitarity requirement similarly as in 
Ref. [23]. Including these virtual corrections and the single scattering contribution, we can rewrite 
the semi-inclusive tensor in terms of a modified fragmentation function Dq^h{zh, tJ^), 



dZH 



= Y,j dxf^{x,t/j)Hl°J{x,p,q,M)DQ^H{zH,fi^) (40) 



where f^{x, 11^) is the quark distribution function which in principle should also include the higher- 
twist contribution [41] of the initial state scattering. The modified effective heavy quark fragmenta- 
tion function is defined as 

Dq^h{zh,H^) = Dq^h{zh,(J''^) 

+ k 4 + (l-.)^M^ ^ X„ -M,-^,,{WlAt.M )Dq^h{zh/z) 

+ i il+%M^ ^ 'fl^^,^,,{z,x,x^A,M^)D,^H{zu/z), (41) 
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where Dq^h{zh, l-t^) and Dg^nizH, IJ'^) are the leading-twist fragmentation functions of a heavy 
quark in vacuum. The modified sphtting functions are given as 



A'yg^gg{Z,X,XLjT,M ) 



/ 



[4 + (l-z)2M2]3^jA(2,^^2) 

A-fg^gg{z,x,XL,iT,M'^) = ^Iq^ggCi-- z,x,XL,eT,M'^), (43) 

Ar//(x,4,M2) ^ I'dz-^ [2Tg^g^Q{x,XL,M^)U=, - {1 + z')T^''gQ{x,XL,M^)] . (44) 
Jo i- — z 

The above equations are similar to the case of double scattering of a light quark [23,26] except that 
the splitting functions for gluon radiation of heavy quark are quite different from the one in the light 
quark case. 

To numerically calculate the modified heavy quark fragmentation function and study the energy 
loss of a heavy quark we have to estimate the twist-four parton matrices T^g'-^{x,XL,M^), which 
are in principle not calculable and can only be measured independently in experiments similarly as 
parton distribution functions. Nevertheless, with some hypotheses they can be factorized. Here we 
will apply the approximation adopted in Refs. [22,23,39], 

^ dy^dy^e'-^P^y-+'-^P^^y^-y^\A\i^g{0)^F,+ {y,)F+^y^)i,giy-^^^^ - yf) 

^^f^ix,)x,f^ix2), (45) 

Xa 

fjly-fly-Qixip+y~+ix2P'^{y^-y2)±i(xL+(l-z)xM/z)p+y~ 

27r 

x(^|Vi,(0)l^F,+ (y2-)F+"(2/r)V'<,(y-)l^)^(-2/2")%- - Vi) 

xa' ^ ' " 

where xa = l/rriNRA, .fq{x) is the quark distribution inside a nucleus, ,fg{x) is the gluon distri- 
bution inside a nucleon and C is assumed to be a constant, reflecting the strength of two-parton 
correlation inside a nucleus. 

In soft radiation approximation, 2; — > 1, the parton matrix elements T^'^^^\x,Xl,M'^) from 
left and right cut diagrams are suppressed and thus can be neglected in our following numerical 
calculation. According to Eqs. (45) and (46), we have 

Xa 

+ f^i^ - (1 - z)xm/z){xl +xt + {1- z)xMlz)fg{xL + XT + (1 - z)xm/z)\ 

^ c,(z,4,M2)C f [f^{x + x:.)xTf^{xT) 

Ixa 

+ fq{x - (1 - z)xm/z){xl + Xt + {1- z)XM/z)fg{XL + Xt + {I - z)xm/z)\ 

- 2fg{x - (1 - z)xm/z){xl +xt + {1- z)xM/z)fg{xL +xt + {1- z)xm/z)] 
H /„ (a; - (1 - z)xM z) 

XA ^ 

X {xL + XT + {1- z)xM/z)f^{xL + XT + {I - z)xm / z) . (47) 
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To further simplify the calculation, we assume xm{^ — z)/z,xl <C xt <C x. The modified parton 
matrix elements can be approximated by 

T^j,-^{x,XL,M^) 

« — f^{x) [(l-e-(^^+(i-")^"/")'/"^)ai(z,4,M2) + a2(z,4,M2)l , (48) 

where C = 2CxTf^{xT) is a coefficient which should in principle depends on and xt- Here we 
will simply take it as a constant. In order to simplify notations we have defined 

ai(z,4,M2) = ci(z,4,m2)-c2(z,4,M2) , (49) 
a2{z, 4, M^) = csiz, 4, M^)/2 . (50) 

Under these approximations, the only parameter in our calculation is C which is also the only 
parameter that enters into the modified fragmentation functions for light quarks [23]. In the study 
of experimental data on modified light quark fragmentation in DIS off nuclear targets, a value of this 
parameter C ~ 0.0060 was found [16] to describe the data very well. Such a value is also consistent 
with that extracted from the study of transverse momentum broadening of Drell-Yan processes in 
p + A collisions [42] . In our following numerical calculations of nuclear modification of heavy quark 
fragmentation functions we will take the same value. 

According to Eq. (41), both heavy quark and gluon fragmentation functions contribute to the 
modified heavy quark fragmentation function. In order to include the contribution of gluon frag- 
mentation, we have to consider the fragmentation function for heavy quarks in the next-to-leading 
order pQCD calculation. Here we follow the ansatz in Ref. [35] and express the overall fragmentation 
function of a parton i into the hadron H as [34,35]: 

Di^H{zH,IJ) f -Df{z,,l)DQ^H{ZH/z), (51) 
Jz„ Z 

where Df{z, fj) is the perturbative fragmentation function (PFF) for a massless parton to fragment 
into a massive heavy quark Q within pQCD cascade. The perturbative fragmentation function (PFF) 
satisfies the normal DGLAP QCD evolution equations [38] and Dq^h{zh / z) is a non-perturbative 
fragmentation function, describing the transition from the heavy quark to the heavy meson, e.g., 
the non-perturbative charm quark fragmentation function into D meson in Eq. (55). 

^Prom the next-to-leading order pQCD calculations [31], we can extract the initial conditions of 
PFF's for the heavy quark at a scale /xq of the order of the heavy quark mass M as 



d:^{z,i,o) = S{1-z) + 



F 



2n 



1 + ^' /^o 



1-z \ 



log^-21og(l-z)-l 



(52) 



Df{z,f^o) = ^ii|l^[.2 + (1 - .)^]log^ , (53) 
^,^,,q(^'Mo)=0. (54) 

With these initial conditions for PFF's and the DGLAP evolution equations we can obtain the 
PFF functions evolved up to any scale /i > /Iq. After convoluting PFF with the non-perturbative 
fragmentation function Dq^h{zh / z) in Eq. (51) we can get the fragmentation function Di^H{zH, IJ) 
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in vacuum, which can be apphed to calculate the heavy quark production cross section within QCD 
factorization formula such as Eq. (6) [35,40]. 

The non-perturbative fragmentation function for the charm quark to fragment into D meson in 
vacuum can be parameterized in Peterson type functional form [30] as 

N 



(55) 



z[l-z-i-e,/(l-z)] ' 

where N normalizes Dc{z) to / dzDc{z) — 1. The parameter Sc is related to the heavy quark mass 
{Mc for charm quark) by — A^/Af^ and A stands for a hadronic scale. 

The PFF at scale /i > /io can be given by solving the DGLAP equations, which is complicated 
in numerical calculations and currently there is no parametrization forms available as for the light 
quark fragmentation function [36]. To simplify the numerical calculations, we choose /io Q as the 
first step for numerical calculations and then obtain D^{z^ Q) and D^{z, Q) according to Eq. (51). 
This approximation is similar to the Approximation Case B used in Ref . [37] . Shown in Fig. 2 as the 
solid line is the charm quark fragmentation function into D meson at = 10 GeV^ after including 
higher order pQCD corrections. These vacuum heavy quark fragmentation functions will be used 
as input in our numerical computation of modified charm quark fragmentation function in Eq. (41). 
We note that in the limit of vanishing quark masses the massless parton model expression should 
be recovered in principle. 
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FIG. 2. Charm quark fragmentation function into D meson in vacuum (solid line) and inside a nucleus 
(dashed line). For the nuclear modification, xb = 0.08 and xa = 0.05 are used. 



Shown in Fig. 2 as the dashed line is the modified fragmentation function of a charm quark into 
the D meson. The value xa — 0.05 corresponds to a nucleus with a radius Ra = 4.25 fm. We 
have taken the charm quark mass M — 1.5 GeV. One can see that the modification due to the 
double scattering in a nucleus for heavy quarks is quite different from light quarks [16,23]. This is 
mainly caused by the form of heavy quark fragmentation functions in vacuum which peak at large 
z. Because of the multiple scattering and induced gluon radiation, the position of the peak of the 
modified fragmentation function is effectively shifted to a smaller value of z. As a consequence, the 
heavy quark fragmentation function remains unchanged, or even slightly enhanced for a large range 
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of fractional momentum z as shown in Fig. (3) by the ratio of modified fragmentation function to the 
vacuum fragmentation function. The modification only becomes significant and the fragmentation 
function is suppressed at large z above the position of the peak. This is in sharp contrast to the case 
of modified light quark fragmentation functions which are suppressed relative to the vacuum form in 
a very large range of z. Note that the heavy quark fragmentation function is strongly enhanced at 
very small z, similarly to the case of light quark fragmentation, due to heavy quark pair production 
from the radiated gluons induced by multiple scattering inside nuclei. However, this enhancement 
is limited to much smaller z than for light quarks fragmenting into light hadrons. 




0.2 0.4 0.6 0.8 1 

z 

FIG. 3. Modification factor for the charm quark fragmentation function in a nucleus. 



V. HEAVY QUARK ENERGY LOSS 



Another advantage of studying multiple scattering of heavy quark in medium is that one can 
actually measure the heavy quark energy loss by flavor tagging, since the leading quark will remain 
the heavy flavor which is very unlikely to be absorbed by the medium. Similarly to the study of 
light quark energy loss, we define the heavy quark energy loss as the energy fraction carried by the 
induced gluon. 



T 



Substituting the approximate expression for the nuclear twist- four parton matrix in Eq. (48), we 
obtain. 



(A.f)(a=B,M^) = ^^ Cdz 
J-^cXa Jo Jo 



[4 + (1 - Z)2M2]4 

(l_e-^'M)ai(^,4,M2)+a2(z,4,M2)] . (57) 

Here we choose the factorization scale as = and define xl = xl + (1 — z)xm/z. Note that 
the virtual correction in Ajq^qq does not contribute to the energy loss. Also, xl/xa = L^/vf 
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with = Rattin/p'^ the nuclear size in the chosen frame. The second term proportional to 
a2 corresponds to a finite contribution in the factorization limit. This term will survive in the 
limit of complete LPM cancellation when double scattering acts like a single scattering for induced 
gluon radiation. We have neglected such a term in the study of light quark propagation since it is 
proportional to Ra, as compared to the R\ dependence from the first term due to the LPM effect. 
In this study we have to keep the second term for heavy quark propagation since the first term 
with the LPM interference effect will have a similar nuclear dependence when the heavy quark mass 
reduces the gluon formation time for low energy heavy quarks. However, for energetic heavy quark, 
the mass can become negligible and one should reach the limit of a light quark energy loss. 
To elucidate the two different limits, we examine the phase factor in Eq. (57), 



{XL + (1 - z)xM/zf _ xl \i\ + (1 - zfM-'f x\M^ 



(58) 



x\ x\z^{\ - zYQ^ ' 

and define it as T = a;^M^/a;^(5^, which should control the LPM interference efifect, and therefore 
the behavior of the total heavy quark energy loss. There are two distinct limiting behaviors of the 
energy loss for dififerent values of of a;s, and xa- 
When T » 1 for < or xb > xa, we have 

~2 / 2 

1 _ (.-^lI^a ^ 1 ^ 

which means there is no LPM interference [20] and we obtain 

Since xa = ^/itlnRa, the heavy quark energy loss in this case depends linearly on the nuclear size 
Ra as the Bethe-Heitler form in Abelian gauge interaction. Similar results were also derived in the 
generalized opacity expansion method [19]. 

In the opposite limit when 3> or <C xa, the quark mass becomes negligible and T 
will take a moderate value. The gluon formation time can be much larger than the nuclear size and 
therefore the LPM interference effect will dominate again. In this case, one can make a variable 
change xl — > xl with xm = (1 — z)xm/z and = fj,^ /2p'^q~z{l — z) + xm in Eq. (57) and obtain 

2, CCAalxB [\ l + Z^ r",^ (XL - XM? 

X [(l-e-^'/^^)ai(z,4,M2) + a2(z,4,M2)| . (60) 

This form is very similar to the one for the light quark energy loss [26]. Since ai{z.C^,M'^) and 
02(2, i^fi, M^) arc dimcnsionlcss coefficients, and T have a moderate value, the exponential factor in 
the above equation coming from the LPM interference regularizes the integration over and limits 
XL < XA- We obtain j dxL/x\ ~ '^/xa- With similar argument in Ref. [23] we can conclude that 
the heavy quark energy loss is proportional to 
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Therefore, the heavy quark energy loss has a quadratic dependence on the nuclear size when 3> 
(X xq <^ xj^. This is understandable because under such condition, the heavy quark becomes 
again relativistic and should behavior like a light quark in terms of induced energy loss. 

To present our numerical calculations of heavy quark energy loss, we rescale the heavy quark 
energy loss by C{Q'^)CAa1{Q^)/Nc and define 

d = ( Azf ) ^ ^ . (62) 

The Ra dependence of the rescaled heavy quark energy loss is shown in Figs. 4-8, where the points 
are numerical results, the dashed lines are linear fit to the numerical results and solid curves are 
quadratic fit. They clearly show that the heavy (charm) quark energy loss has different nuclear size 
dependencies as and xb change. In Fig. 4, Fig. 5 and Fig. 6 we fix — 10 GeV^ and change the 
value of Xb to demonstrate the nuclear size Ra dependence of heavy quark energy loss for a charm 
quark (Af = 1.5 GeV). One notes that when a;^ is very small (large heavy quark energy), the charm 
quark energy loss depends quadratically on Ra (see Fig. 4). However, as we increase xb (decrease 
heavy quark energy) a gradual transition from a quadratic dependence on Ra of the energy loss to 
a linear dependence (see Fig. 5) takes place. When is very large the charm quark energy loss has 
an all most linear dependence on nuclear size Ra as Fig. 6 illustrates. Similarly, shown in Figs. 6, 7 
and 8, we see that even for large values oi xb (where we fix xb = 0.15) or small heavy quark energy, 
the same transition from linear nuclear size dependence to quadratic dependence of the heavy quark 
energy loss takes place as we increase . 
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FIG. 4. The Ra dependence of heavy quark energy loss for a charm quark. 
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FIG. 5. The Ra dependence of heavy quark energy loss for a charm quark. 
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FIG. 6. The -Ra dependence of heavy quark energy loss for a charm quark. 
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FIG. 7. The Ra dependence of heavy quark energy loss for a charm quark. 
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FIG. 8. The Ra dependence of heavy quark energy loss for a charm quark. 

Another mass effect on the induced gluon radiation is the "dead-cone" phenomenon [18] that 
suppresses the small angle gluon radiation. Since the size of the dead-cone 0o = ^ll" [Eq. (34)], 
within which the gluon radiation is suppressed, is inversely proportional to the quark's energy, the 
reduction of energy loss is stronger for a slower heavy quark. For a heavy quark with either a high 
energy or virtuality Q^, its radiative energy loss should approach that of a light quark. Since 
the dimcnsionlcss coefficients a\{z,i\,M'^) and a2{z,i\^]VP') in Eq. (60) also depend on the heavy 
quark mass, they will have addition mass effects on the heavy quark energy loss. 

To illustrate the difference of energy loss between heavy quark and light quark and the quark 
mass effect we define a ratio R as: 

_ (AzQ)(2;b,^2) 
^= (Az|)(xB,/i^) ' ^^^^ 

where (A2;^)(xb, A*^) is the light quark energy loss [26] which can be obtained by setting M = in 
Eq. (57) and Eq. (60), 

(Az,«)(.i.,M^) = -^ r ^ rdz(l + z^)(l-e-i/-i)C^[l-i-i]. (64) 

In Fig. 9, we show the change of the ratio R with for a charm quark (Af — 1.5 GeV) propagating 
in a nucleus with xa — 0.04 and xb = 0.1. Please note that in the nuclear parton matrix elements, 
the fractional momentum in a nucleon is limited to < 1 due to the momentum conservation. Even 
though the Fermi motion eff'ect in a nucleus can allow > 1, the parton distribution in this region 
is still significant suppressed. Thus it provides a natural cut-off for x^ in the numerical integration 
over z and It in Eqs. (57) and (64). 
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FIG. 9. The dependence of R for a charm quark. 

We can observe that when is not too large, the heavy quark mass effect significantly suppresses 
the energy loss caused by induced gluon radiation. When A'P /Q"^ 0, the efli'ect of quark mass 
becomes negligible and R ^ 1. This is consistent with the pQCD factorization theorem that when 
the momentum transfer is very large one can neglect the effect of quark mass. 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 
FIG. 10. The xb dependence of R for a charm quark. 

Shown in Fig. 10 is the Xb (or heavy quark energy) dependence of the ratio between heavy quark 
and light quark energy loss for fixed = 10 GeV^. It is obvious that suppression of heavy quark 
energy loss due to the "dead-cone" effect of heavy quark mass is most significant when is large 
(or quark energy is small). When is very small (quark energy is large), the effect of quark mass 
is small and the quark energy loss approaches that of a light quark. 

VI. SUMMARY 



Utilizing the generalized factorization of twist-four processes we have studied the nuclear modi- 
fication of heavy quark fragmentation functions and the energy loss of a heavy quark propagating 
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through dense matter after it is produced via a hard process in DIS in the twist expansion approach. 
Taking into account of the multiple scattering suffered by the heavy quark we have derived the modi- 
fied heavy quark fragmentation functions related to twist-four corrections with nuclear enhancement. 
We find that the formation time for gluons radiated from a heavy quark is smaller relative to that of 
a light quark, since it is always measured against the propagation time of the quark. With certain 
kinematics when the quark energy and virtuality is small, the gluon formation time can become 
much smaller than the nuclear size. In this case, the heavy quark energy loss or the nuclear modifi- 
cation of the heavy quark fragmentation functions have a linear nuclear size dependence. We have 
shown through both analytic deduction and numerical calculation that a gradual transition from a 
linear to a quadratic nuclear size dependence takes place when one increases the quark's energy or 
initial virtuality. 

We also compared the energy loss of a heavy quark with that of a light quark and demonstrate that 
the quark mass effect, including the so-called "dead-cone" phenomenon, will significantly suppress 
the heavy quark energy loss when the momentum transfer is not too large. This heavy quark mass 
effect will decrease if the heavy quark energy, or the momentum scale is much larger than the 
quark mass. When M — > 0, our calculations recover the results for massless quarks in previous 
studies. 

Similar to the case of light quark propagation [16], the results discussed in this work can be 
easily extended to a hot and dense medium, which will have practical consequences for heavy quark 
production and suppression in heavy ion collisions. When the data on direct measurement of D- 
meson spectra in high-energy A + A collisions become available in the near future, one should be 
able to use the modified fragmentation function in a parton model to study the modification of 
the D-meson spectra [43] and probe medium properties similarly as one has done for high px light 
hadrons [15]. The different pattern of energy loss for heavy quarks, such as energy and medium size 
dependence, will not only confirm the unique feature of non-Abelian energy loss but also give more 
confidence in using jet tomography to study properties of dense matter in heavy-ion collisions. 
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APPENDIX A 

In this appendix we list the calculation results of double scattering of the heavy quark discussed 
in Section III. There are total 23 cut diagrams as illustrated in Fig. 11-21. 
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There are three three different cuts(central, left, right) in Fig. 11, and their contributions are 



Hi {y ,yi ,y2 ,kT,x,p,q,M,; 



, (l + z^)4 + (l-z)^M^ 

^ (1-Z)(£2, + (1-^)2^2)2 

X -^Cf ^^Iifiiv^ ,Vi ,V2 J'T,kT,x,p,q,M,z) , (Al) 

X I _ p-i(xL+{l-z)xM /z)p'^y:^ ^ _ ^-i{xL + {l-z)xM /z)p'^ {y^ -v:[) (A2) 

/i,L(2/",yr,2/2",^T,A:T,:r,p,g,M,z) = -e^'^^+^-^^^y' -y^^ e{y^ - y^)e{y- ~ y^) 

liAy- ^vi .tTM:X,p,q,M,z) ^ ~e^^-+^-^^^y-+^--^^^y^-y^^e{~y^)9{y^ 
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FIG. 12. 



In Fig. 12, there are two different cuts for induced gluon radiation, central or left. They give, 



-TpD 



H2{y ,yi ,y2 ^kT,x,p,q,M,z) 



(1 + Z^)^T ■ (^T - (1 - z)kT) + (1 - zfAP 
(1 - + (1 - zflVP][{fT - (1 - z)kTf + (1 - z)2M2] 



^ 27r ' Y'^^ ^^ '^1 '^2 jT,kT,x,p,q,M,z) , 



(A5) 



/2,c(y-,2;r,y2^^T,fcT,x,p,g,M,z) = e^(-+--)^'+^"+^--^'+(«i"-^^")^?(-y2-)%--yr) 
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-i(x L + {l-z)x M / z)p+ {y -y-^ ) _ -i{x ^ + (1- z)x m / z)p'^ (y -j/i +2/3 ) 



(A6) 



-i{xL + {l-z)xM/z)p+{y^-y^)+i(x'j-,-xo)p^{yi ~y2) _ Q-i{xL + {l-z)xM /z)p+ {y^ -y^ ) 

where x^j = /c|n/2p+g~. 
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FIG. 13. 



As for the central cut and right cut of Fig. 13, we obtain 
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FIG. 14. 
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There is only one cut (left cut) in Fig. 14 with the contribution, 
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FIG. 15. 



As for Fig. 15, we get 

Hfiy^,yi,y2,kT,x,p,q,M,z)^J dej. 
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FIG. 16. 



The contribution from Fig. 16 is 



Hgniy ,yi ,y2 ,kT,x,p,q,M,z) / dl 



(1 + Z2)4 + (1 - ^)4m2 
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X QiixD—-^D)p'^{yi-y2)g^-'i'i^L+{l-z)xM/z)p+y^ _ (A16) 




FIG. 17. 

As for the processes in Fig. 17 we have three possible cuts with the corresponding contributions, 



Hr^ciy ,yi ,y2 ^kT,x,p,q,M,z) = / dt. 



(l-z)[(^T-fcT)2 + (l-z)2M2]2 



^ ^^^~^^7,c(y ,2/1 ,^2 jT,kT,x,p,q,M,z) , 

ZTT I\c 

77My',yi,y2,f^T.kT,x,p,q,M,z) = e^^-+^-^p^y'+^--p^<^y^-y^^e{-y,)e{y- 



(A17) 



X e 



ia:Dp+(j/~-j/j")/(l-z) _ -i(a;i, + (l-z)2:M/z)p"'"(j/"-aj") 



(A18) 



H^MR)^y ^Vi ^y2 ^kT,x,p,q,M,z) = J 



^(l-z)[£2, + (l-z)2M2]2 



X — ^?^/7,L(fl)(y ,2/1 i^T,kT,x,p,q,M,z) , 



27r A^c 

-i(l-2:/(l-2))a;£5p+(j/j"-a2") 
]^ „ g-i{xL + {i-z)xM /z)p'^ {y^ -y^) 

l7,R{y~>yT>y2,eT,kT,x,p,q,M,z) = -e^^-+^-'>p^y-+^^-p^^y^-y^'>e{-y^)9{y^ ~ y^) 

y. g-i{i-z/{l-z))xDP^iyi -y2) I _ g-i{xL + {l-z)xM/z)p+y2 



(A19) 



X e 

X 



(A20) 
(A21) 
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FIG. 18. 



There are also three cuts in Fig. 18 and we have 



(1 + z^)£t ■ {It - kr) + (1 - z^M^ 
(1 - + (1 - z)2M2][(£t - kr)^ + (1 - z^A'P] 



as Ca 27ras-= , _ _ _ . 
2tt~2 >yi ,y2 jT,kT,x,p,q,M,z) , 

78,c{y-,yi,y2,eT,kT,x,p,q,M,z) = -e''-^+--^'p^y-+"-^^^y^-y'^^e{-y^)e{y- ^ y^) 

X I _ Q~iixL + {^-z)xM/z)p+{v^ -y^) 

78,L{y',yi,y2.^T,kT,X,p,q,M,z) = -e^(-+--)f+«"+"-''^fei"-^2-)0(j^- 

-i{l~z/{l-z))xDP'^{yi -y2) 

^~i(xL + {l'~z)xM / z)p^ (y^ -y^) _ ^ixoP* (y^ -y^) / (l-z) 

78My',yi,y2,eT,kT,x,p,q,M,z) = -e^(^+^-^p^y-+^^-p^^y^-y2)ei~y^)9{y^ -y^) 

X Q-^(^i' + ('^-'^)^M/z)p+y2 _ ^ixDP^y^ /0--Z) 



X e 

X 



(A22) 



(A23) 

(A24) 
(A25) 




FIG. 19. 



The contributions from the three cuts in Fig. 19 are 



-ttD 



Hg{y ,yi ,y2 ,kT,x,p,q,M,z) = / dl 



(1 + z^)£t ■ {It - kr) + (1 - zfNP 
' (1 - z)[ej, + (1 - zYAP][{£t - kTf + (1 - z)2M2] 



otg Ca 27ras— , _ _ _ 

— ~J^^^{y '2/1 '2^2 ,iT,kT,x,p,q,M,z) , 



(A26) 
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-i(xL + {l-z)xM / z)p+ (y^ -y^) _ ^ixoP^ [y^ -y^) / i^-z) 

79My-,yi,y2JT,kT,x,p,q,M,z) = -e'(-+--)f+^"+"°p'-(^r-^=")0(-y2-)%2- -yf) 

X g-''(l-z/(l-^))a;r)p+(Kj"-i/2") 

X p-i(xL + {i-z)xM/z)p'^y2 _ ^ix DP* y2 / i'^- z) 



(A27) 

, (A28) 
(A29) 
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FIG. 20. 

^From Fig. 20 with two possible cuts (central or left), we have 



Hio,c(.y ,yi ,y2 ,kT,x,p,q,M, 



Ml 



(1 + z'^){iT - kr) ■ [tr " (1 " z)kT) + (1 - zfA-P 
(1 - z)[{fT - krf + (1 - zfM^][{fT - (1 - z)kTY + (1 - zYAP] 

as Ca ^Tvasy ^ 
X— — -j^iw,c(y ,yi ,y2 ,iT,kT,x,p,q,M,z} , 

ho,ciy-,yi,y2JT, kT, x,p, q, M, z)e^i-+--yp^v-+^-opHy^-v^) 9{-y^)e{y- - yf) 

y.g~i{xL + {l-z)xM /z)p'^y2 ^ixDP'^iy^ -y^)/{l-z) _ ^-~i{xL + (l-z)xM /z)p+ {y^ -y^) 



(A30) 
(A31) 



7?^ 



w,Riy ,yi ,y2 ,kT,x,p,q,M,z) 



(1 + z'^)iT ■ (It - zkr) + (1 - zfAP 



(1 - z)[tl + (1 - z)2M2][(^T ~ zkrf + (1 - z)2Af2] 
x-T^ ^^^^Iw.R{y^,yi,y2jT,kT,x,p,q,M,z), (A32) 

^-i(2;r>-a:n)p+(a7-y2^)-«(2:t + (l-z)2;Af /2)p+yi~ 

^-i(l-z/(l-2;))2;Dp+(a7-a2")-*(2:i.+(l-^)3:M/2)p+t/2^ ^^33^ 
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FIG. 21. 



Finally, Fig. 21 has both central cut and the left cut with contributions as 



Hiiniy ,yi ,y2 ,kT,x,p,q,M,z) 



(1 + z^)iiT - kr) ■ [It - (1 - z)kT) + (1 - z)^NP 
(1 - z)[{fT - krY + (1 - zYNP][{fT - (1 - z)kTY + (1 - zfM-^] 

as Ca 27ras— , _ _ _ 

r^^ii,c(2/ ,yi ,2/2 ,I^T,kT,x,p,q,M,z) , 



(A34) 
(A35) 



Hii.Liv iVi ,y2 ikT,x,p,q,M,z) 



(1 + z^Yt ■ {It - zkr) + (1 - zfM'^ 



(1 - + (1 - z)2M2][(£3, - zkrY + (1 - z)2M2] 
X ^ ^^^ii,L {y^ ,yi,y2,£T,kT,x,p,q,M,z), (A36) 

ZTT Z l\c 

~hi,L [y- , Z/r ,y2,iT,kT,x,p, q, M, z)e^(-+^-)P^y-+^-oP^ iy^ -vl ) eiy- - y-)9{y^ - y- ) 

^-i{xD-x%)p+{y^ -y^)-i{xL + {l-z)xM /z)p''' {y^ -y2) 

-i{l-z/{l-z))xDP'^{y^ -y2)-'i'{xL + {l-z)xM /z)p''' {y^ -y^) (A37) 
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